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VORTEX F FORMATION AND ND RESISTANCE 1 IN PERIODIC s 

"McNown! M. ASCE and G. H. ‘Keulegan 
d 
; a In two recent, independent studies, (1, 2) closely related aspects of the 
esistance | of a bluff body to the periodic motion of the ‘surrounding fluid were 


 neailangeel. In combination, the two studies provide a clear distinction be- P. 
_ tween effects due to velocity and those due to acceleration as well as some in- 
sight into their variations during unsteady flow. Although indirect, the ae 
of viscosity is paramount in that its consequences are separation, vortex 
formation and shedding, and resultant alterations of the virtual mass. The 


_ specification of these various aspects provides a basis for the correlation of of 


_ The effects on resistance of the variation in size of a wake during unatenny >. 


has already been discussed in some detail.(2) A motion which starts a 

from rest is initially irrotational and is free from separation. As the velocity _ 
increases, edge vortices form and grow into a wake. The completeness . 
the wake | formation depends upon the comparative ‘magnitudes of the du duration — 
of the unsteady flow and of the interval of time required for a wake to; form, 7 
and upon the degree to which the body is streamlined. il a ia 

If f the form of the wake | varies with time » the c corresponding flow patterns © 

and pressure distributions around the body vary. The forces attributable to. 
pore or to acceleration must then also change with time. . Because both _ 
kinds ¢ of forces depend upon the state o of development of the wake, they are in 
turn related to each other. For example, while the wake of a bluff body  . 7 
small, the large curvatures of its boundary results in low pressures onthe 

_ +rear of the body and in a comparatively high coefficient of drag; simultaneous - 


ly, ‘the disturbance to the flow is less and the coefficient of mass is compara- _ ™ ei 


Fe Experiments at both the National Bureau of Standards and at the University 
of Michigan \ were planned to disclose the c characteristics of resistance in un- 
Steady flow. standing wave wave was produced in an open tank so as to 
Note: oc open until June 1, 1959. To extend the closing date one month, a “— d 
___- written request must be filed with the Executive Secretary, ASCE. Paper 1894 ee ; 
f part of the copyrighted Journal of the Engineering Mechanics Division, Proceedings _ 
tk the American Society of Civil Engineers, Vol. 85, No. EM 1, January, 1959. | 
a. Presented at the October 1957 ASCE Convention in New York, N. Y. 
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two- -dimensional oscillatory motion past a cylindrical in 


_ nodal plane of the wave. Thus, the flow varied from zero to a maximum a 
one . direction, back to zero again, and then reversed direction, so that simple — 

_ harmonic motion was closely approached. By varying the amplitude of the 
wave and the size of the body, one could obtain a ‘motion with a comparatively | 
“small amplitude and negligible separation at one extreme and quasi- steady 
tad at the other. The resulting occurrences can be distinguished by means a 


ofa dimensionless frequency parameter UmTo/D, in which Um represents 


; transverse dimension of the cylindrical body (measured normal to the direc-_ 


The work at the National Bureau of Standards was based on the “eer YW 
of motion (see List of Symbols), 
T p dU 


a [kau U ak 


in two ‘respects. combines the effects | ate pressure 
7 gradient and of inertia. Also, the effects of the term U dk/dt are absorbed, 
for Eq. (1), in the variation of Cm. The two equations differ conceptually in 


gration was feasible in either case, the inner boundary for the contour inte- 
es gral was the surface of the body for Eq. (1), whereas the e boundary wasa 
7 * surface ¢ enclosing both the body and its time- dependent wake for Eq. (2). sea” 7 
a two studies. In the first, an average value was obtained for a given occur- 4 
_ Fence by harmonic analysis of an entire cycle. * In contrast, at Michigan the 
_ hypothesis was utilized that the value of k varied with time during each | “u 
cycle because of its | dependence upon the size of the wake. Accordingly, : 
: combination of the two methods of approach | indicates only general trends in 
which average values are compared with instantaneous ones. 
‘The free-streamline theory of Riabouchinski for flow around a flat plate 


_ Was extended to provide a theoretical relationship between 1 drag and virtual = 


~ mass.(2) This resultant was used with an assumed variation of the virtual 
_ mass | coefficient with time to calculate the e cyclic variation of to total force. sh. 
i - Although | evidently assumptive, the approach provided strong evidence of the i 4 


importance of sizable and related variations of the coefficients of drag and 


Results from experiment for flows past a thin flat plate and a circular 


In the two parts of Fig. 1, Cm is | 


| plotted a against Cp. ; The | points for the cate with the flat plate and the aver-— 


Barameter from Figs. 10- 13 of the NBS study. The plot is striking evidence 
of the existence of unique relationships between the two parameters. Also, 
the theoretical result for t the flat at plate is a | a surprisingly good first t approxima- 
_ Because time has been eliminated as a parameter in 1 Fig. 1, the plots are : 
& one ‘step removed from e physical o occurrences. halt alf cycle at 
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_PERIODIC MOTION 
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36 
“ Inter -relationship Between Coefficients of Drag and of Virtual | - 


ass for (a) Flat Plates and (b) Circular Cylinders. 
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moment of zero ambient to some point at the upper 
end of the extended curve in Fig. 1a. _ The ensuing occurrence is indicated by 4 
a partial traverse of the curve, the terminal point depending upon the relative 


duration of the half cycle, and a return along the same curve to the ‘starting / 


the second half of the cycle. The numbers indicate relative time in a dimen- | ¢ 
Sionless form described in a paragraph. Because each experi- 


during each cycle was undoubtedly even greater than it appears to be in Fig. i. 
a Both parts of the figure indicate the existence of marked ‘secondary | ef- 
fects which are ascribable to the formation and shedding of vortices andto 
= the existence of vortices created in the preceding half-cycle. The interpreting 
of the graphical indication is somewhat presumptuous because the necessary — 
concepts are still in the formative s stage, and because the methods of ai analysis 
- Establishment of flow past a flat plate is greatly affected by the sharp _ 
edges ‘which lead to separation for very small relative ‘motions. “Initially, a 
_ symmetrical pair of vortices forms and grows at the edges, so that regions of 
_ ‘negative pressure also form : and d grow. In ac accordance with the predicted ri re-_ 
sult, the smaller the vortices the lower the pressure along the downstream 7 


fee the plate. At first, these do not extend across the entire downstream 


and one of the vortices begins to o grow more rapidly than the other, , eventually ; 
becoming | so large that it breaks away from the plate » and moves downstream. — 
During the period of growth, the vortex is still attached to the body, and Cy 
: increases markedly as S the composite body becomes fatter. ‘Then, as the vor-— 
‘toni is shed, _ the net form is again more sylph-like and the effective value of _ 
‘Ce decreases abruptly. Both trends are readily apparent a ‘Fig. la as the 
a two. marked departures fr from the predicted monotonic curve. ‘The ‘shedding of 
other vortices produces comparable effects, but, ‘because the data were 
gy ponevtow- for each cycle, only the effect of the first is clearly delineated. ais 4a 
motion in one direction persists long enough, in a relative sense (i. if 
fi: -CmT/D is large), a succession of vortices is shed as for a taut wire singing 7 
in the wind—the well-known phenomenon of the Karman vortex trail. The! flow 
is then quasi - -steady. The average value of Cp decreases t toward the ac- 


“boundary decreases, and the value of Cy increases with increasing 


average § size of the wake to a value of about 5. 


a parity. _ The replacement of the sharpness 0 of the edges of the flat plate “2 
= the comparatively generous rounding of the circular cylinder eliminates 
= = the initial separation and the consequent regions of very low pres- 
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the limiting value of unity tor Cm for some large but finite value | 
ie 
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q 
™ jj. an alternating vortex trail is established. This orderly progression is well a 
documented by the photographi us portions of 
in Fig. 1b arre not asr 
— 


available. The minimum value occurs, as before, as first vortex 
shed. Also, because the vortices do not extend the disturbed region 
laterally, as they did for the plate, the value of Qi does not increase dur- 
- the period of initial growth. The lower branch ¢ of the curve for the cylin- 
der is comparable to the lower portion of the curve for the plate. = = 
The extent to which vortices develop depends upon the duration: of flow in 7 
“one direction. ‘% More precisely, the ratio of this time interval to a time inter- 
val characteristic of vortex formation is a fundamental parameter. Use of a 
_ the Strouhal number to indicate the time to form one vortex pair in a com- = 
B =, steady flow is advantageous even for - unsteady flow. Thus, the oot 
- o of the motion is some multiple (or fraction) of the period Ty of vortex o at 
= maximum velocity). If is or r less, separation and vortex 
tion are relatively unimportant, the inertial effects being approximately a 
for the classical unseparated flow, and the | velocity | so small that t the drag 
force is negligible. If the ratio is 10 or more, the motion becomes quasi- 
on The numbers around t the curves in Fig. 1 are values of To/Ty. In each La 
~ case the e complete period of motion is used (a full period for the wave, and 
" the time for a pair of vortices to be shed). For several reasons, the num- 
< bers are 2 or 3 times as large as they should be to fit exactly the pattern —_ 
_ described. The effective velocity is obviously less than the maximum, and 
the maximum occurs at the quarter-cycle. Vorticity from the preceding — 4 
cycle noticeably retards the formation of the first vortex. Also, the insta-_ 
; bility which causes one vortex to be favored over the other of a pair may be 


flow. _ Whatever the cause, a study of motion pictures of the phenomenon con 
‘firmed the supposition that vortex formation is delayed for the periodic mo- Ps 
tion studied, and that the primary characteristics of the curve are directly _ 
related to the growth and shedding of the vortex as ‘Suggested | in the preceding 


paragraphs in spite of the | numerical disparity. 


‘slower to develop in flow which starts from rest than it is in steady snc 


7 _ The foregoing concepts of resistance in unsteady y motion are offered as 


a 
basis for further investigations. Its several shortcomings and cffered as a : 
can be remedied only by the completion of additional studies. Significant E 
Progress should result from the evaluating of the effects of growing and mov-— 
7 ing vortices. These are predictable from classical theory, and can be cor a - 
related with measured forces. . Studies should be devised to test and to ‘extend 4 
_ these concepts so as to interpret rather than merely to describe the effects 
of unsteadiness. So many kinds of unsteady flow are possible that a wholesale 
“empirical solution is impracticable. . The separation of the effects of velocity vi 
and of acceleration, the designation of the relative time scale appropriate to 
phenomenon, and the definition of the of vo are essen- 


force} force per unit length of cylindrical 


"a ‘inertial coefficient (including term for p pressure gradient) — 
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Oscillating Fluid,” Natl. Bur. Std. Rep. S. department 


J.8., “Drag in in ‘Unsteady Flow,” IX Int. Appl. Mech., 
Prandtl, L., and Tietjens, O. G., Applied Hydro- | and Aeromechanics, 
York, McGraw- -Hill, 1934. 
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APPLICATIO NS OF ELECTRICAL ANALOGS OF STATIC STRUCTURES 


certain similarities between elastic in structures and power in electri- 
Bi networks, and in which current is used to simulate forces and moments — ; 
in the structure, and voltage to simulate ne and p rs, 


"continued mere in following 


2. Treatment of structures containing members which are 


either by reason c of a non-linear stress-strain curve for the material 
or because the deflections are large enough to significantly alter the ‘ 

equilibrium relationships among the forces and moments. 
3. A practical ‘method of ‘setting up analog ‘circuits. 
4. Partition of complicated structures to permit simulation by an ana a 
device with a limited | number of electrical 


- first appear, in| the illustrations or in the text, , and are arranged ‘alphabetical- 
Bi ly, for convenience of reference, in the Appendix. Except for new symbols wa 


oe herein, the notation is the same as in the companion paper.\"’ (1) 


‘Wow: Discussion open 1, 1959. To extend tt the closing date one 
2 ait _ written request must be filed with the Executive Secretary, ASCE. Paper 1895 z= 
Ae part of the copyrighted Journal of the Engineering Mechanics Division, Proceedings 
i the American Society 0 of Civil Engineers, Vol. 85, No. EM 1, January, 1959. abel 
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The basic "approach | followed by the writer in electrical analogs 


for static structures has been explained in the Introduction of the companion > 
paper. (1) , Since the writing of coat paper other works on the subject have ap- 
peared. (2) _In general, however, are b based on the consider- 


_ Hence these papers do not deal with structural energy and electrical pov power — 
explicitly, except possibly to make use of the former in evaluating the de- ae 


flection and slope values needed in step 2 above. By contrast, in this s and the | 
_ companion paper(1) the second and third of the above considerations are com-_ 


The « ‘consequence of requirements 1 4 can 1 be realized 
by comparing Castigliano’s well-known theorems with their electrical equiva- - ‘s 
 lents(3) which have been described in the companion paper, as follows: a 


Castigliano’s Theorem Electrical Equivalents 


Forces and moments are as Currents are such as to = 
to minimize the elastic energy, electrical power, subject to re- 


subject to. ‘quirements of current- continuity. 


The deflection (or slope) at : a given 2. ; The voltage at a given point is the 
energy with respect toacorre- 


sponding force (or moment) ‘applied ‘Supplied to this this point. 


— 


work in the structural case) the forces and moments take on definite Shag 


(often designated simply as Castigliano’s Theorem in the structural case) the 


, and the currents take on corresponding values. As a result of the second ~ 


Ps deflections and slopes at all —_— take on definite values, and the voltages _ 


_take on corresponding values. Thus by measuring the voltages in the 


4 _ flections, slopes and voltages were not considered explicitly (in steps 2 and 3 


The companion paper deals with the for a number of illustrative 


_ types of structures. The present paper treats various topics which were Jett 
out of the previous paper in order to keep it at an introductory level, _ -: 
Bp tetsu the material had not been developed at the time that paper was on . 
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STATIC STRUCTURES 
Dimensional Frame 


oA treatment of the for the 3- dimensional 
is too lengthy for inclusion here but may be found in the thesis (4) on which 
this | paper is partly based. Instead we shall investigate a single specific 
which however displays most of the features of interest. ie ae 
‘Fig. 1a) shows three orthographic projections of a rigidly connected frame 
in which the cross-section is constant in any one member, and which is load- — 
ed by force F and completely restrained with respect to moment (including . 
_ torsion) and ‘force at support points 1 and 4. Fig. -1(b) shows an oblique pro- - 7 
jection of the frame and defines the coordinate axes x, y, and z, and also 
_ Shows the three unknown forces and three unknown moments at each of the two 
- support points. The moments are indicated by circled arrows, each of which E> 
a : denotes the direction of advance of a right-hand screw rotated by the moment. 
_ There are a total of twelve unknown reactions and six equations of equilibrium 
p of the overall structure (three force equations and three moment oqustions), — 
that it is indeterminate to the sixth n degree. 
In developing the analog we shall make use of the moment analog 
_ shown in Fig. 4(c) of the companion paper,(1) and reproduced here as Fig. 2b) 
and in abbreviated ‘symbolic form as Fig. 2(c). | The corresponding structural 
_ member is shown in Fig. 2(a), where for our = 
sidered that the shear S and the respective applied end-moments MBH and a 5 | 
Myup lie in the plane of the paper. It is seen by inspection that the oon 
of moment equilibrium of the beam is simulated by current-continuity re- ae 
quirements in the analog, and it has been shown in connection with Eq. (15) in 
the companion paper(1) that the electrical power dissipation (with the trans- — 


former assumed ideal) will equal twice the it 


where E is the modulus of elasticity, 1 is at 
the beam cross- section, and the expression on the right is the bending com- 
ff _ pliance, or change of slope per unit constant moment in the beam. (It is as- 

_ sumed that the work associated with direct shear and with any tension which — 
may exist in the beam is negligible.) It has also been shown, as a consequence 
of Castigliano’s Theorem and its electrical equivalent, that the voltage drop 

from B B to H in Fig. 2(b) represents the increase . of slope between B and H- 


‘ae deflection of H from the line tangential to the beam at B, inthe di- 
rection o of aw shear; and correspondingly for the voltage drop from H to : 
The suales of Fig. 1 is shown in incomplete form in Fig. 3. It is cided 7 
‘into three ‘more or less separate networks, as indicated by the horizontal dot- 7 


lines, dealing essentially (but not exclusively) with X-, y-, and Z- 


in Fig. 1(b). ‘First considering member 1- -2, » we ground point 1 in all three 
aa to simulate zero values of x x-, y-; and z- slope,* and designate the : 
values of current leaving the networks at these points to ‘correspond with the — 


_ Slope is here defined as the elastic rotation in the direction “F a given com- 


me of moment, regardless of whether the rotation arises from flex 
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amount Z Lo between 1 and 2: this is taken care of by using the analog meee: q 
in Fig. with a suitable value of shear current. Corresponding consider- 
= hold for z- moment in member 1-2 acts in torsion rather than bending, 


‘ where G is the shear modulus and I7y is the moment of inertia, effective in 
computing torsional deflection, of the beam cross-section. The voltage drop I 
across this resistor in the direction from 1 to 2 is then the increase of y- < @ 

» 


slope, which corresponds with the fact that the voltage drops from 1 to 2 in 
- x- and z-networks refer to x- and z-slope respectively, = 
aa In using the” bending analog for member 1-2 we have assumed that the | 
cross-section of the member is symmetrical about at least one of the two | » 
axes to which the components of bending moment have been referred, in this — 
_ Case the x and y axes, else the resulting slope changes could not necessarily — 
have been determined separately; a corresponding assumption will be made 
for the other members. Of course the values of resistors within the moment | 
analogs corresponding to the two components of bending moment for the — 
_ member will depend on the values of I referred to the corresponding axes. y 
ae By inspection the currents designated as Mx23, My23 and Mz23 equal the “ae 
--—s« corresponding components of moment in the structure (as a matter of oo- 
nition, Mx23 is the component of moment, applied to end 2 of member 2- 3. re-. 
ferred to the x- -axis). The three electrical analogs for member 2-3 are hag 
eae as before; the x- ~moment produces torsion and the zZ and z- -moments 


[In dealing with member 3- 4, ‘whose length is Wt. it will be convenient to 

_ refer the moments to axes x, u and v v, where the latter two are designated by 
arrows in Fig. “For y and z are also 


_ Moment Mx34 produces inten. and the aie increase in mome 
- between 3 and 4 is (Y, + Z4) Lo; this is taken care of in Fig. 3 by a moment 7 
analog as before. My3q is a torsional moment, , and is handled (within the y- 


_ network for convenience) by a single resistor as before; the means of develop- 


_ a corresponding increase of - y 2(F + Xj)Lp along the member, and is handled 
accordingly. Points 4 are grounded in = to simulate zero slope 


ing the current M,34 will be dealt with later. My34 is a bending moment with 4 q 


corresponding unknown reaction moments defined in Fig. 1(b). Moment Mx; 
electrical power dissipate 
| 
— 
a ¢ 4 _ spectively. The moments Mx34, My34 and Mz34 then represent correspon img 3 7 
components of moment applied to member 3-4 at 3, and the voltage of point 3 
a 


ry 


3 


(incomplete) 


cm — 
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Ou3 | and Oy3. if desired, the y- and 3 
termined as follows, as is evident by inspection of 1a): 


@73= Gy3)/V2 
on The 1e complete analog is shown in Fig. 4, which is obtained from Fig. ‘3b 


formers and external current-generators as necessary. primary and 
a of any one transformer are designated by circles connected by a 
1) ) and the required turns ratio (includ- 
ing polarity) is fixed by the ratio of primary and secondary currents. Any 
ni scheme of connection which satisfies the designated currents is acceptable _ 
4 with this proviso: : the various currents simulating any one unknown force a 
7 moment must be so inter-related, either by transformers or by direct con- 
nection, that if one of these currents is determined the others are fixed by | - 
a ‘considerations of current -continuity ; alone (here extended to include the effect - 
of the fixed ratio between the primary and secondary currents of each trans- 
former). For example, current -Z\Lo at point Q is related to the current — f 
"Balla at P P by direct connection, and to the corresponding current at R by . fi : 
ad transformer. _ If this transformer were omitted, the number of redundant ol 
_ currents (those not determined by considerations of current-continuity alone) 
would be greater by one than the number of redundants in the structure, and — Pe 
‘the analogy would be invalid, 
Zz 7S Castigliano’ s Theorem and its electrical equivalent, the voltages at Q 
and O, , measured below ground, are respectively inn times the mepeisroe “a 
629 and -4x2 of joint 2 in the z and -x directions (the y-deflection is sero). 
the voltages at Pp and are 6. v3 


“librium: in the structure. The summation of are leaving the x- ‘portion of 


the y-portion two current summations may be 


e currents leaving the circui "in the y- and z-networks, 
the cur _ The voltages at points 3' in the y- and z 
— 
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STRUCTURES 
wherens in the z-portion the following ng two may be ead 


pe My y34 a" Mz34 may be eliminated by Eqs. (6) and (8), and My4 and My4 a J 
mi may bs expressed in terms of My4 and Mz4 by inspection of Fig. 1. When the 
resulting substitutions are incorporated in Eqs. (7) and (9), andthese 


4 equations are successively added together and then eel one from the 


| 


_ which are seen to be the muinien of equilibrium of y- 


_ Experimenta work in connection with this example is discussed elsewhere | 


Non- “Linear Structures 
— linearity in structures may be divided into two general classes: the — 
— one arises from a non-linear relationship between stress and | strain, and ie. 
includes the case where the stress- strain relationship is not unique but de- — 
| Struc on the previous loading history. Important applications of this occur in 


elastic limit but before complete destruction has taken place, as in . bomb 
_ ghelters. In the second class the stress- -strain relationship may be linear, 


conjunction with a large axial load to and 


In a pin- \-connected truss this may be treated by preparing for each member a 
a curve of compliance vs. axial force, and then by successive iteration find- — 
ing such values of those resistors which ‘simulate the compliance (as in ‘Fig. .- 
of the companion paper(1)) as are compatible with the axial force values. — aa 
a practical matter the compliance is known in advance for relatively light 
loads, so that a convenient method of procedure | starts with the analog set for 
light load, the latter then being increased in small steps so that the successive i 
_ adjustments in the resistors are small and hence relatively independent of 
each other during any one step of the iteration. The method may be extended ns, 
to vary the load (in small steps) according to any desired loading program, _ 
watch is useful when the compliances depend on the loading history. ae 
In a containing bending the situation is much more 
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_ lates =~ and the vertically- disposed resistor may be adjusted so that the 
voltage drop from B to O simulates 6pH/L. As in the previous | case, it is 
well to start the analog ‘at light load and to change the load in small steps. 3 


Non- Due to Effect of Deflections on Equilibrium 


This o occurs in individual members it in rigid frames. and 
ses, and in practically all members of long slender trusses subjected to 
“forces which tend to place the structure as a whole in compression, as in | 
= boom trusses. In individual column-like members it is often possi- EC 
ble to neglect the extra deflections and slopes due to column effect until the — 
load is large enough to cause column failure. On the other hand, when the en aa 
tire truss acts as a column the effect of the deflections on the force- “we ar 
equilibrium equations may have to be considered. As a highly simplified ex- 
_ ample consider the pin-connected truss of Fig. 5(a), where it is assumed that — 
the structure. as a whole is s adequately restrained against column buckling 2 
_ outside of its plane, and that individual compressive members have no ap- . 
preciable column instability. Conditions of symmetry permit us to split the 
truss along the longitudinal axis of member 3- -4, and to consider the left i 


the changes of slope of the various members caused by the relatively small 
x- and y-components of deflection taken into account, the transformer turns 


ratios must be such as to impose suitable ratios between the x- and y- 


components of tensile force inamy member, — 
__ The procedure is to first set the transformer turns ratios with the de- == 
flections assumed to be zero, find the resulting deflections by measurement 
of the analog, and then compute new values < the turns ratios. Successive 
runs are taken until the solution converges, or (if the truss is unstable) di- = - 
4 verges. _ However, the transformer voltage drops can be expected to be rela- 
tively very large since they have been shown to represent x- and y- | wie 


- components of branch deflection due to branch rotation, which hare ew 


imperfections in greater detail below) will be 
large, and the analog isn not recommended for this type of application. -<elgl 


— 


Practical Method of Setting up Analog Circuits 
the case of complicated structures it is convenient to have a relatively 
a =a practical cases is described in the following Feng in connection with a 
_ the 2-dimensional rigid frame of Fig. 6, where it is assumed for saapiictty Ca 
that cross-sections are constant, and that the elastic energy associated with | - 


change of slope along the member in Fig. 2(aj, as well as the deflection 
at H measured normal to the tangent through B are functions(5) of both 
a Mpy and For any values of these two parameters the horizontally- 
| 
- 
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| at 4. The x- an -networks oi the electrical ana T Fig. ar 
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5 COMMN-LIKE TRU TRUSS AND ANALOG 
be used, as well as nore a : 
Choose reference directions for axes x and y as well as for slope (9), 
and show the loads and reactions, the latter assumed unknown. For ae 
uniformity show reaction moments in the direction of positive slope, 
resolve the reaction in the x- and where 


a 3. Assign ¢ a reference positive | direction, and t 


alphabetical direction of the joint designations; if the opposite direction 
eg is preferred as a reference in some cases, indicate this by an arrow on 48 
- - Designate the x- and y-forces acting on each member, as dictated by dig 
equilibrium considerations alone. Use the following conventions: above 
orto to the left of the member indicate the x- -force applied to the member a 


4 

= 

| 


STRUCTURES 


e 


pan 


6 RIGID FRAME 


"sponding y- -force . The result is shown in Fig. 6. Where 
7 forces are redundant in that they cannot be obtained by equilibrium con- 
_ Siderations alone, express tl them | in terms of the unknown tension 


when it tends to rotate the member in the direction of positive 
-_ a oy For example, fo for a horizontal member such as d-g running ~ 
from left to. right, a force at the beginning of the member in the positive 
aia ore -direction is the negative of shear force.) The joint-force equilibrium» 
considerations must also satisfy the two force summations 1 al 


cow 


atYe+ 


q 

— 
4 
y 


‘calculations since it permits one to think of the x- or y- 
electrical or hydraulic flows in the positive | direction of each member. 
Draw a new picture of the frame p) preparatory to. developing the analog 
Ye in Fig. 7(a). Ground the terminals where the slope is held to zero, 

and leave open the terminals where the moment is zero. Drawthe _ 
analogs of individual members, making use of the symbolism of - oP 
Fig. 2(c) wherever "(The analog for member e-h may be 


out affecting total power, to take advantage of the fact that one ne end is a 
moment-free). Compute the appropriate shear currents by inspection | 
of the x- and y-forces acting at the beginning of each member, and Hj oe wre 


_ designate these values but leave the shear-current ‘connections open. ; 
6 . Now use ingenuity freely in setting up transformers, current- => he 
generators, and inter-connections SO as | to develop the ‘necessary 
currents, ; as shown in Fig. Tb). The three equations of overall equi- 
Librium must also be satisfied. (10) is satisfied by the currents. 
approaching point Z; Eq. . (11) be used to develop acurrent Yo 
(which does not appear in the analog), but this is not done here because _ 
it is so easy to perform the equivalent operation by direct computation — 
be after the other unknowns are determined; and the equation ofoverall 
aon ‘moment equilibrium is satisfied automatically, : as can in be seen by | taking g 
the summation of all currents leaving the 
which the moment summation about joint d 
‘There are many other ways in which the shear currents can be de- | ia 
‘a veloped. All are satisfactory if the conditions of redundancy are identi-— 
cal in the analog and in the structure, as discussed in connection with | 


7 tf 7. Assume a unit of length which will give a convenient value for some © - 
length the structure. Choose convenient scale 


aaa for each ‘member, “switch on the generators, and read nee 
_ Finally, to check the results of the work by an entirely independent ~~ 
method, determine by ‘structural analysis’ whether the values of force 
and moment yielded by the analog are in equilibrium, and whether they 
result in slopes and deflections compatible with each other. In general 
there will be one compatibility check for each degree of redundancy. 


the present case these equations are satisfied by the x-forces and y- 
: 
current to moment, The scale factors between voltage and slope, and 
j- between voltage and the ratio of deflection and L,, will then be fixed. 
4 
weer 4 q 
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: In the practical use of the analog, cases will arise where the sonden « of 
available electrical elements (such as transformers, resistors, and current 
- sources) is not sufficient to permit simulation of the complete structure. we 
: _ means of simplifying the structure so as to overcome this limitation will be 
7 described in connection on with the particular case of ot Fig. 6, which has been — 
‘discussed above. ‘The n method applies only | to linear structures, in which the : 
response to successive equal of a given force or moment 
remains the same. 
First we cut the structure into two -structures of which at 
least one is in static equilibrium (though not necessarily statically determi- 
ate as shown in Fig. 8(a). At the inter-connections of this portion (in this" 
case the right-hand sub- -structure) with the e remainder r of the structure there 


are unknown forces and ‘moments as shown, the symbol F being used to indi- sy 


cate generalized force, which can be either true force or moment, deocies 
the following equations | giving the generalized deflection 5", y which denotes ra 

translational deflection when associated with true force, and rotational ae ; 
flection or slope when associated with moment, the prime indicates that intl 


where F5 and F6 are ignored since they are not siansiinhed with deflection. | 
= constants a ajj are conventionally known as influence coefficients. By | 


where W' is the energy of the ~hand -structure e with Q' ignored. 
trial we find that this is compatible with previous expressions if 


whens use has been made of the reciprocity theorem, which states that 


i 
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(a) subd- structures 


= 
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. Continuing this technique in successive steps we finally 

_ By inspection we determine that it is in general possible to chenae values _ 
for the forces such that all but any one of the square terms equal zero, 0 Since : 


each square term must be positive, SO that it can be simulated by an electri-_ 


ing power dissipation will be proportional to twice the elastic. energy of the ca 
-sub-structure. If in addition these resistors are connected to an analog which 

_ simulates the left-hand sub-structure in Fig. 8(a) in such a manner as to ote 

‘simulate equilibrium requirements by current-c ontinuity considerations, the 4 

_ analogy will be complete insofar as this portion of the en is concerned, ; 

ot 
To take the latter into account we must withdraw from the sila an 

; amount of power equal to twice the work, W", which is imposed on the right- 7 

»@ hand sub-structure by the displacement of t the generalized inter- connection = 
- forces through the generalized deflections 6" caused by Q' (or in general, by 
7 all external applied loads acting on this sub- structure). To find 7. imagine 


their final values, in which case the inter-connection forces and their associ- — 


= all external loads on the structure are increased uniformly irom zero to 7 2 | 


deflections also from zero to their final values, sO 


| 


sources. To find the values of 6" iv as well as of the ‘constants in Eq. Ga) a a 
consequently Eq. ( 

Fig. 8(a) is an actual ‘structure, set up its electrical pio determine 
the inter-connection deflections with the various inter-connection forces and 
external loads (in this case Q') successively applied. 
Asa preliminary to drawing the analog, we denote x- and y-forces acting — 
ni at the beginning of various members of the left-hand sub-structure in the > 
Same manner as was done in Fig. (6), and then express the inter- -connection 
forces in terms of — forces and moments om considering equilibrium at at 


where the denend on certain of the a’s_and Fy anneare ie 
— 
4 ealactrical rent flowino t ouch the eonrreenondinge resictor then the ecult- 7 
7 
A 
4 
— “al 
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ee: As to considerations of overall equilibrium, moment equilibrium will pre-_ 
sumably be automatically : satisfied by the analog as before, while force = 
‘item yields no additional information over - what we already have. jis 7 
_ The analog is shown in Fig. 8(b), where the portion representing the left- -s 
hand sub-structure is developed by the same method as b before. _ The currents. 
Fga and F4 | appear directly in this portion of the analog, , while F3Lo. and Fil, = 
5 are developed with the aid of current generators Gj and Gg in accordance 
with the above equations. Power equal to 2W' is simulated by means of eteene- 
formers and resistors in accordance with Eq. (16), where however the coef- 
ficients of F4 have been factored out wherever F4 appears between parenthe- 
ses; this saves several transformers by avoiding the need for multiplying > 1 i 
. by various values. The power represented by Eq. . (17) is simulated with the 
aid of voltage generators of suitable polarity, as indicated, to absorb ‘power 
since W W" represents work absorbed by the right-hand sub- 


represents summation of the left- hand sub-structure ure about 
In the present case it is ‘clear that these is no great saving in circuit com- 
ponents resulting from the use of the partitioning technique. Obviously this _ 
-method increases in usefulness as the number of members in the first su oem 
increases | and the number of inter- -connection forces both 


agpliied aad ts, structures are often subjected 
n addition to applied forces an moments, structures are often subjecte 
: to applied deflections and slopes, as when a joint of the structure is displaced © 
or rotated by a prescribed amount. This is closely associated with thermal & 
_ effects, which change the (unloaded) shape of a member, so that stresses re- 
sult if the change of shape is restricted by adjacent members. _ This phe- oa 
4 nomenon is important in piping systems carrying fluid at varying ‘tempera- 
oe general r method for modifying the analog so as to deal with applied chang-— 
es of shape, whether resulting from thermal effects or 1 not, is given in the | 
_ following rules. Like the rules relating to the remainder of the analog, these 
employ current- ~continuity and energy concepts alone, and their validity - 


easily proved by the energy theorems previously ‘stated, me the reasoning» 
: _ 1. Select a force (or moment) F' which would result in the deflection (or 
_ Slope) 6' of interest, and arrange the analog so as to simulate F" by a 
aa 2. In series with this current insert a voltage generator, of voltage | pro- - a 
ee - portional to 6', the constant of proportionality being such that the ower 


supplied to the analog by the generator equals F'6'. 
4 pve an init suppose that joint 1 in Fig. 1 is given a known deflection am: 


in the y-direction, or alternatively that member 1-2 alone is heated to such e 


a temperature that its free expansion would be a. The structural effect is — 


the same in either case, the associated force being Yj. To simulate this it is - 
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necessary to at point S in Fig. 4, where the current Lo flows, 
a voltage difference 6'/Lo, with polarity such that the voltage rises in the _ 
direction so that the power Yj 5' is supplied to the circuit. All 
currents and voltages in the analog then automatically adjust themselves aa 


_A new model of the equipment has been built and is shown in Fig. 9. ete 
need not be described in detail because it is essentially the same as the Pu ; 
equipment shown in Fig. 11 of the companion paper. (1) However, the — 


arrangement. 


othe incorporates a 400- -cycle power which 
_ smaller errors due to transformer imperfections than the 60-cycle _ 
supply previously used. ¥ Further, a quadrature- rejection circuit of a type 
conventionally known to electrical engineers is incorporated in series with a 
- the panel meter, which, as before, acts asa naps scp in combination with 


eel . is the sainaioa’ effect of the series and shunt inductances of the transformers, 
The series resistances of the transformer windings can often be absorbed 
sao within analog resistances to which they happen to be connected. Sometimes, Ss 
_ however, ‘such contiguous resistances do not e exist in the analog, st so that it may 
es be « desirable to compensate for the winding resistances, as well as for the ef- 
fective shunt resistance of each transformer, which in general cannot be ab- 
_ sorbed within analog resistances. This can be effected with the aid of power 
admitted to the transformer through : a compensating winding, which may or a 
may not be part of the primary or secondary; however, no such compensation 


= 


= 


Where more than one transformer carries a a 
single force or moment, their combined winding resistances may be referred 
to a a single winding in n such a way t that the total power is unaffected. 5 ll 


i In the. following experimental w wok each of the force errors in a given 


_ structure is is expressed not a as a percentage of of its own theoretical value, but me ie 


= asa a percentage ofa single reference force selected as the largest a 
— most significant in the structure. The reason for this choice of error cri- 


terion is that, because of practical imperfections both in the analog and in the : 


‘ 


actual structure, the various forces in the analog and in the structure tend to 


_ differ from their theoretical values by absolute values which are of the same _ 


|" 2 
— 
4 
— 
4 _ ation than in the previous case, and space has been conserved and accuracy 
"a a _ increased by omitting the resistor knobs and dials and instead providing for ; im 
g of the resistors using a built-in Wheatstone 
24 
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order of magnitude throughout the structure, aan which are not highly de- 
pendent on the magnitude of individual forces. Thus, errors referred to indi- 


vidual values might ergen very large in the case of lightly loaded members, | 


and hence are easily taken into account by the designer through a a ‘slight abso- Ps oa 
i ——_ lute increase in the strength of lightly loaded members over and above what - 
is theoretically required. This type of increase would be needed even with ee ey; 
perfect analytical solution, because of geometric imperfections and un- 
certainties ofloading ofthe structure. 
i A _ corresponding argument holds: for moment errors. Also, by similar 


reasoning, deflection and slope errors may be expected to tend toward the £ 
same order of magnitude throughout the structure. | 
'* To avoid confusion due to signs, an error is called positive or negative ac- 

cording as the absolute magnitude of the associated experimental value is - 

larger — than that of the theoretical value. 


_ _ Members 1-2 and 2-3 were assumed to have a circular cross- ss-section, so a 

tha 


t the flexural moment of inertia I is constant for all components of bending 
moment, and the torsional compliance is 5/4 times the bending compliance 
(with the shear modulus assumed equal to 2/5 of Young’s modulus). "For 
-member 3-4 the values of I were assumed | the same, but the torsional com- - 
_ pliance per unit length was considered twice as great as for the others. " i 
Fig. 4 was modified by multiplying the voltages and dividing the currents Ps 
_ in the sub-network ABCD by v2, thus permitting current X Lo to flow direct- ; 
ly from point O to point E, E, and current My34 from points F and G to point H. ae 
r Transformer shunt resistances were ignored, and winding 1 resistances were 
"absorbed within analog resistances except for transformers handling 
34: and was taken as unity, Lo/EI as 100 ohms, and V2 Fa 


Frame of Figs. 4 = 


-0. 
0.3347 


on 


il 
| 
q 


All cross-sections are , assumed equal, and load Q' is is twice as large as Fo 


or Fe, which are equal. Transformer shunt resistances were ignored and 
winding resistances were absorbed within analog resistances. Lo was a 
- as unity, L,/EI as 100 ohms, and Fp as about .025 ampere. In the following — 
table the theoretical values for forces are factors of es for moments of ely 


rable Frame of Figs. 6, 7 
Theoretical Analog Theoretical Analog 7 


Unknown. “Value _Error Unknown — Error 


(71.6248 
-0.3752 
-0.3739 
70. 
0406 
0.1900 
0.1900 


Partitioned Frame of Figs. 6, rig 


_ To obtain a good test for the validity of Fig. { 8(b) it was considered sila 
able to determine the values of the influence coefficients in Eqs. (14) and (16) 
_ by analysis, rather than by the application of the analog to the right-hand sub- 
structure in Fig. 8(a) as would be done in practice. Table Ill shows pertinent — 
values including those of use in constructing Fig. 8(b). With L, assumed — 
2 equal to unity, the numbers given for influence coefficients are » factors of 


/EI. The table also shows values of 6", with the numbers representing 


‘The factors between ‘compliance and resistance and between 
’ and moment were the same as in the previous case. Transformer shunt re- 
‘sistances v were ignored, and id winding resistances were compensated either by © 
being absorbed within analog r resistances or by suitable adjustment of the 
_ voltage sources simulating 6". Since the resistors in sub-network ABC were 
a relatively very low, the voltage and currents in this sub-network were -— ae 


spectively multiplied and divided by 10 by y means of a suitable connected 
‘Theoretical values and analog | errors relating to certain unknowns, inc 

_ ing the generalized inter-connection forces Fj through F4, are given neat 
“Table Ill, where the factors of appropriate quantities as a 


‘Frame of Figs. 

d to X, and XaL,, respectively, and the deflec 
moment errors are referred to X, 
4 
— 
Ez : 
3 
— 

“4 


Values Needed for Setting Up Analog 


208 boy 


Value ___Err Error, % gar ‘Value Error 


is usually to derive the governing is 
- quired with general-purpose analog computers), and it is particularly easy to 7 
explore the effects of modifying the dimensions of structural members, , which 
e: is useful in design procedures. The use of energy theorems make it possible a : 
to obtain slopes and deflections without considering them explicitly in setting — 
; a up the analog circuitry, resulting in in relative simplicity a and in economy of 
a” The relative ease of veiuiitie: a separate ‘tai makes it feasible tode- _ 
tect gross errors due to misuse or defects of the analog. 


pen to such things as conditions of end fixity, dimensions of members, joint © a 
conditions, and the neglect of direct shear and tensile effects in bending woe. 
‘members. . Thee error may become significant in special cases, such as in — 
shallow trusses subjected to high longitudinal compressive load, where the q 
voltage drops are very large the voltage drops 
Techniques have been included for dealing with applied deflections and 
en effects, and for partitioning structures which have too many membe - 


simultaneous simulation on an analog instrument. 


— 
q 
i 
Z 
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influence coefficient; 
cross-sectional area; 


F force associated with deflection; current analogous to 
same; 


shear modulus 


moment of inertia. effective 


- moment of inertia. effective in torsion; 


load; current to distributed load; 


fo force; current nt analogous to tensile 
"rectangular coordinates; are 


"generalized ¢ deflection ca by external loads; 


elastic slope; voltage eatemen to e 
resistance analogous to compliance. 


1. “Blectrical Analogs of Statically Structures”, by F. L. Ryder, 


2. "See for example “The Direct t Analogy Analog Computer,” by G. D. McCann, 
Journal, April” 1956 which | an extensive 


Transactions ASCE, Vol. 119, 1954, Pp. 1046. 


Network Analysis by Least Power Theorems, 


“The 


i 
- 
— 
ae: 


A . treatment of this is is given in “Strength of Materials, ‘Part 
Timoshenko, D. D. Van Nostrand C Co. New York, 1941, p. 371. 


6. ‘See to Higher 


_ New York, 1907, p. 131. 

2. Similar results be obtained by dealing with in a manner 

scribed by Cauer in “Ideal Transformers and Linear Transformations, 

Elektrische Nachrichten Technik, Vol. 9, p. 154, May 1932; however, the 7 
el present method is more appropriate here in that it is based on energy con-— 

ie in 1 place se of the compatibility considerations used by Cauer, as : 
a consequence of which the transformer voltage drops need not be con- © 
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ABSTRACT 
a 
of longitudinal waves in an elastic infinite plate 


_ is discussed by applying the Method of Internal Constraints and by taking in 
-account second order terms in the equations of constraints. Numerical ~<a 
_ Sults obtained by applying this theory are compared with those obtained by 
applying the exact theory given by Lamb. we 


LATURE 


NOMENCLATURE, 


The following ng nomenclature is used in the text: wat 


v,w =C components s of the elastic displacement in the x, x, y, zaxes 

of strain 


thi ‘thickness | of the plate 


ara. 
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4 : a. The results: presented in the paper were obtained in the course of an © 


investigation sponsored by the Office of Ordnance and by the 
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INTRODUCTION 
_. The problem of two-dimensional waves in a solid bounded by iin = 
planes has been considered by Lord Rayleigh(1) and by Professor H. Lamb(2) 
in two well-known papers ublished respectively in the years 1888 and 1917, a 
- Recently, J. P. Marsden( (3) numerically computed dispersion cu curves for longi- 
_ tudinal waves from the theory given by Lamb. The object of this paper is to — 
"discuss the problem of sin papal longitudinal v waves in an infinite ==. 


m- 


a one- -dimensional approximate theory. of wave propagation based on 

_ the Method of Internal Constraints. It was the late Professor R. M. Davies — 
= advised the authors to compare the ‘results given by the approximate oF ; 


e authors are indebted to the late Professor R. M. Davies, former Head ‘'¢ 

q of the Department of Physics at the University College of Wales in 
Aberystwyth, for having furnished them with numerical data which have been 

used for comparison in some of the figures enclosed in the text. gues 


— 

— 
— 

— 

_ The present paper complements two previous papers,‘\*,°/ in which the 

7 _ problem of dispersion of longitudinal waves in elastic straight rods of a a 

— 

— 

field vectors and their variations over the eross-section.(6) 

Pa _ Ca The infinite plate of constant thickness h will be referred to a system of ‘a 


Oo | 


0 


Calling u, v, w the ‘components of the in the x, ® 


2) normal nal 


E being modulus and s ratio tio oft the of the bar. 


Eq. (3) into the ‘components of the elastic dis- 


On the a assumptions that no no extereal forces (body an and surface forces) ar 
applied to the infinite plate and that the plate has constant thickness h, the ow : oe 


— 
hat on the free surfaces o! ing equati—_ 
= __ The condition that on | s expressed by the following as wists eee 
, Should vanish is exp | (2) 
— 
— 
> 
— 
eo | 
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K which appears in Eqs. (4 Eqs. (4) is defined 


i 
Where A, Ag, 2, Ag are constants, phase velocity, the wa 
‘number, L the wave and i = 
; 
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By substituting Eq Eqs. (6) into Eqs. (4), the following | three homogeneous 

age 


ae 


By expressing that the determinant of poe 


3; 


(4), the ratios cand ‘in determined 2 as follows: 


— = — 
+ | | 
| aa 
| 
a 


v value of 


be used. . Expressing t the factor kK by Ec ‘Eq. i in of an, one 


According to the Exact Theory (Lomb) 
= K, 


to the Present Method, kK 


to the Method, K= = _K, 


| 
= 


3 

| in order to determine the 
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"LONGITUDINAL WAVES 3 = 


>0.29. This been plotted in Fig. 2 and has compared 
the exact theory given by Lamb(2) and computed by Marsden. (3) 
4 ‘This root becomes complex for a value of c close > to the value of t the veloci- 7 
ty of the Rayleigh surface waves. Apparently, the reasons for this are: : first, 
= the parabolic distortion of the cross section assumed by Eqs. (1) is not valid — 
for ‘Short wave lengths, and second, the influence of the shear stress on ‘the | 


paper, does not vanish. 

e ratios =~ an 


“this ca case Eq. (12), ac closer approximate value for the factor K has been ob- 


have been later computed for K - 21, in 


a? 


root of Be (8) for SS has be 


ie eaneiapanl corresponding numerical values so obtained are plotted in in Fig. 2. The ¥ 
agreement with the exact theory given by Lamb seems to be excellent. zi Also 
. 
in this case, however, this root becomes complex for values of c close to ad ‘iz 
values of the velocity of Rayleigh surface waves. _ 
_ Components of Elastic Displacement and Stresses Along = Cross — of 


for the components. of the « elastic seplacemen’ 


which has been designated as Ko. The variation of with is shown 


= 


+a ) (se ¢ 


ig 
j 
| 
| 
ca 
Wwe 
P 
= 


0.95 90 85 0.8 0.75 


Correction factor K asa function on of 


In | (13) and (14) A= iA). 


— 

— incorre- 
ore. 
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ot In Fig. . 4, the ratio of the longitudinal displacement at the eaxtece of the 
plate u(h//2) to the ‘displacement at the is plotted as a 
of the ratio Inco rrespondence to the value id nodal 
surfaces appear on the two surfaces of the the piste. ch 
In Fig. 5, the nondimensional _— 


7 Ratio of displacement at the surface of the plate to the saieeaee 


at the middle surface as a function of a 


' 
| 


4 


Variation of displacements and stresses over 
‘the: cross-section ofa — of thickness he 


F I UR RI 
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: Ceci by the method presented : in this paper for a value == 0. 217 using 


are compared with the 1 numerical results computed by J. J.P. Marsden 
fh — - 6, 7, 8, 3, the above ratios are e plotted as functions oe the east 


x 


and for the of particle velocity, the 
in drawing the drawing these diagrams, , the following convection was ii ie the 7 
= _ constant A is positive, the above ratios were taken with the plus sign when | = 
the displacements and stresses near the middle surface oft the were of 


at the form or iA(Re where (Re) denotes a positive q 


real number. if Instead, the » Gaplatoments and iapating were e of | the form 


4 a -A(Re) e and -iA(Re) e the above ratios were ‘elie with the if 
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ariation of the longitudinal displacement = 
over the cross-section of a plate of thickness h. 


— 
— 
— | z= 
= 


a 


that tay and canb be > written under the the — 


= 


~ 


effect of ‘the of the factor K discussed in in 
t ‘the paper. 


a 

| | “| 0.328 | — 

Yt — 
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Variation the transverse stress over the c: 
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_ where V is the volume of the plate. By iaiueleiitine: Eas. (16) and (18) into 

Eq. (20) by performing the integration with respect to the variable y x, 


Where is the length of the in x , the of the 


x By ‘substituting Eqs. (17) and (19), into Eq. (22) by performing the se 


respect to variable the for the kinetic energy 
obtained: 


1-7 


— On the assumption that no external feieee es (body and surface forces) are 
asain to the wae, by applying Hamilton’s Principle under the form: — 


to the whole plate, supposing the displacement to 1 os zero ro at ¥ sy te, the 
equations of motion ( (4) are obtained. 
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ELASTIC STABILITY OF THIN SPHERICAL SHELLS 


Gideon P. R. von Willich 
_ Previous theoretical work relating to elastic buckling of thin in spherical Aon 
; shells is reviewed. A theory, based on strain- -energy considerations, is py 
_ presented for determining buckling pressures of shallow thin spherical a 7 
— _under ¢ external pressure, = it is suggested that the results may be extended — 


During recent years, -shell-type structures have increased i in popularity. 
Where light weight has been sought, the buckling characteristics of these _ 
types of construction have often been an important consideration. It is there- 
os not surprising that the question of « elastic stability of shells has been of | 
great interest to aircraft and missile designers. However, with the tendency a 
in building construction towards the use of shallower and lighter shells, this — eo 

is also becoming important to structural engineers. 
‘Unfortunately, even simple physical phenomena relating to shells tend a 

lead to intricate and tedious mathematical analysis. This is especially the 
case in the investigation of buckling problems. This paper is therefore 
restricted to consideration of a single aspect, namely the elastic stability of f a 


thin under ext external pressure. 


Review of Previous Work 
On page 491 of reference 6, Timoshenko describes the “classical” or linear 
li buckling theory, together with a short history of the problem. In this approach 
Note: Discussion open until June 1, 1959. To Rati the closing date one month, Bsc 2 
pi part of request must be filed with the Executive Secretary, ASCE. Paper 1897 is - 
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part of the copyrighted Journal of the Engineering Mechanics Division, Proceedings _ _ 
the American Society of Civil Engineers, Vol. 85, No. EM 1, January,1959. 
1. Prine. ‘Research Officer, Mechanics Sub-div., National Mech. Eng. _ 
be Research Inst., South African for — and 
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it is re that buckling will occur at a pressure | — permitsan | 
-equilibetem deflected shape that is infinitessimally removed from the trivial 
spherical shape. The method is thus analogous to the determination of the 
Euler buckling load of an initially straight strut under axial compression. ~ a 
_ This method of attack adopts the “eigenvalue” view of buckling, in which it is" 
- supposed that the structure retains its shape, undergoing only uniform con- _ 
traction, until the critical load is reached. ‘ At this point, a deflected shape — 
- differing from the trivial one becomes possible. The linear theory does not © 
; allow the determination of the magnitude of the deflection at any point; with — 
a constant buckling load, arbitrary amplitudes are permissible. In the case : 
of a uniform strut, it is possible to obtain a solution of the governing differ-_ 
ential equation using a more accurate expression for the curvature of the —_ 
strut, and it is then found that the load does increase slowly with increasing 7 
aes displacement, although at the start of buckling the rate of increase is a 
= _ The analogous problem for the shell, with the use of more accurate _ 
= for curvature, is not readily solved, but it can be accepted that 
the equilibrium pressure decreases with increasing deflection, sothat 
the configurations exist at pressures appreciably below the predicted 
_ buckling pressure. In this way, the behaviour of the shell differs essentially 
from that of the strut. Experimental evidence also shows that actual buckling o 


“roughly between two and four times too high. | The discrepancy in n the latter 
case also extends to the buckled shape: the theory predicts that this will be 
in the form of a number of ‘small waves covering the surface of the shell (all 

a of which will be in unstable equilibrium under constant pressure), while in 


"practice buckling is localized to a single “dimple ” which grows progressively 


a: Several papers have been published in attempts to explain the discrepan ae 
te Von Karman and Tsien, (8) proposed the hypothesis that the practical pers 
~ buckling pressure of a spherical shell will correspond to © the 1 minimum pres- 
sure necessary to keep the shell in equilibrium in a buckled shape with finite — 
deflection. _ However, _ Friedrichs(2) showed that certain assumptions made le by ; 


these authors were unjustifiable as they : seriously affected the results « ob- 
_ tained, and that the hypothesis was not tenable. A subsequent ps paper by Tsiet Tsien a7) 


1. energy level (i. e. strain- -energy ir in the shell plus potential | energy y of 


the applied pressure) must be the same before and after buckling, and 


The of the loading procedure must satisfied. 


“tal ‘evidence... _ However, the criteria stated are ‘open to question. energy 


—— there is no reason why it should not be lower. In fact, experimental — ey 7 


observations show that the snap- -through type of buckling in shells is generally x 
| eri by noticeable noise and vibration, indicating a definite energy 


loss. It is of course necessary that the loading conditions must be satisfied, 7 7 


— but it. is difficult to understand why they should affect the magnitude ofthe 
—ueteling load: they merely influence the behaviour of the shell after the maxi- 


load has been reached. s distinguishes" between — 


« € 
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‘pressures correspond closely with theoretical values in the case of struts, s- 
in CASO O herical shells dicted bucklino pressures 
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loaded by conainnt pressure (for instance in structures immersed at a great 
_ depth in a fluid) and those in which the pressure is reduced during testing (as 
is the case when a specimen is tested in a tank with liquid under pressure). ty 
Tests do not show a difference in the buckling presmures obtained | by the 1 two. 
- methods, and one cannot assume that the shell “knows” before it commences =, 
a to buckle, what will happen to the load during the buckling process. gon 
The “eigenvalue” approach was adopted in all these publications. This ap- aK 
proach, with its supposition of discontinuous behaviour at the buckling load, — 
applies to perfect structures, which are not attainable in practice. It is more 
é logical to regard the buckling phenomenon as a loading process in which the | 
: relation between load and deflection is non- -linear. . In the case of the strut it 
*. _is readily shown that the behaviour predicted by the “eigenvalue” procedure | 
[* - is the limiting case as imperfections tend to zero. However, it is doubtful © 
whether this is true for a spherical shell, and the mathematical complexities — 
of non-linear deformations of spherical shells make this question difficult _ 
if not impossible to settle. s 
The ‘differential equations governing axi-symmetrical deflections of thin 
"shells of revolution were given by E. Reissner, (5) the equilibrium equations 


_ and the stress-displacement equations being reduced to two simultaneous non- a 


_ linear differential equations in two variables. The solution of the equations 
presents enormous mathematical difficulties, however. For the case of 
4 shallow spherical shell of uniform thickness t the form of the equations be- 
San much simpler, as some of the terms become negligible. The simpli- S 
fied equations are still difficult to solve, but significant progress has been | 
_ made in obtaining series solutions, for instance by Kaplan and Fung, (3) Reiss, 
_ Greenberg and Keller,(4) and Archer.(1) However, a large number of terms’ 
in the series solutions are required to obtain convergence. = ~~ 
> ‘Several authors have reported experimental results, which are illustrated 
_ graphically in Figure 2. Tests were also carried out by W. Delano in the v 
_ Civil Engineering Department of the Massachusetts Institute of Technology ; 
in 1953-54, under the supervision of Prof. C. H. Norris. re ve 
shallow spherical shells clamped around the circumference, with 
given thickness and radius of curvature, the buckling pressure varies with the if 
_ distance of the circumference from the axis. Ny For a certain value of this dis- - og 
tance, the buckling pressure will be a minimum, and it is felt that this mini- 
q mum is of some significance in predicting the buckling pressure of the gh 


Scope: The following | analysis refers to shallow thin spherical shells; how- 


ever, up to a certain stage, the equations are equally valid for other shallow 
thin shells of revolution, as this extra degree of generality can be conveniently — 

introduced. By “shallow” and “thin” are meant that, respectively, the central a 

; Fise and the shell thickness are small in comparison with the base diameter. =— 


‘The shell is assumed to be rigidly fixed around the the boundary an and to ” sub- — aa 


The thickness is assumed to be constant, and the material to be meee 
a E ous and isotropic, and linearly elastic with the same modulus of elasticity in 
~ tension and cc compression. It is assumed that the elastic limit of the — 
is not apiaaead Deflections of the or are assumed to be symmetrical about 


jected to uniform normal external pressure. 
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ana to 0 the : middle st surface after deformation without extension 


Notation: Some of used are illustrated in the accompanying 
ea 1. Others are defined where they first occur in the text. Cylindrical — 
polar co-ordinates are used. ‘For brevity, di or 


denoted by ac a comma and a a subscript, e.g. 


Strain 


“Strain Expressions 
derived on the assumptions that 


strain at middle surfa surface 


dial in curvature 


“= 


_ ‘The following relations for the three ¢ energy aidan are easily derived: 


train rgy to extension of the middle surface ‘atest 
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= _ Displacements 


Si 


~ 


>: 
hy = circumferential of cur vature 


: 


| 
i 
= _K,= radial change of curvature 
i 
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The princiale of minimum potential can to find the relation be- 
_ tween applied pressure and central deflection. 


TMs principle can be stated: 
Of all displacements satisfying given boundary conditions, t those which — 


% ‘satisfy the equilibrium conditions make the potential energy assume a sta- 7 
= value, and for stable e equilibrium the Potential energy must be a i 
‘ . The procedure adopted here is to to assume a deflected shape, containing 
certain unknown parameters, and to” apply the principle of 


of the mal has two re ay of freedom, , corresponding | to the two — 
ment components. However, in using the strain-energy procedure in connec- 
_ tion with an assumed deflected shape, it is sufficient to assume a Suitable . 
shape for the > component "parallel to the shell axis. — 
_ for the other component @% can then be determined byt 
of the of variations. 

The total | 

U, contains: u, hence given 2/) the condition for minimising 


ait the ante expressions are expressed in terms of displecementa, as ae 


given by equations (1) and (3), the following differential equation is obtained: _ 


For a shell, is a known function of and w is assumed in terms 


| 

rey of the applied pressure is 
4 | 
| 
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— 
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and hence, 


‘If the right hand side of equation (8) can be written as a polynomial in >, _ 


the gene 


a= 


From the nature of th the deflec deflection sian and of z a it is found that 


. 6, in fact that contains only powers. r. 
‘The constants B will be known in any particular | cane, ‘and At be 


Up to this. par the analysis is applicable to thin a shells of revolu a= 
“thon; the subsequent is restricted to segments only. 


shell is to deflect in the form 


‘Tos satisfy the boundary ry conditions at = 


W,,= Q, thisbecomes 


a *+856C+2048) 
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Determination of Critical Pressure 


correct is to minimise wi to the two variables s C 


For Bad : : The equation resulting putting A ow. =0 in 
equation (9) may be written more compactly by setting: aut i 


PAL 
-271 ct, 148.0. + 
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fan) 


e ee value of Ww W being given} by the the smaller Of the | two rv roots. ots. 


ot 


165 1 1 


hon 
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Then 


q minimum value > Of Por “with the corresponding value Cc bys a ti trial and 
For 0, the ‘corresponding equations are as follows: 


= 


2054+ 484074 +5808C* 
+ 25300 


condition for buckling is = 0, which gives = 4 
| 
ting this into equation (11). 
a given value of A and hence of 77 , the optimum value of Ccanbe 
found, in theory, by setting Poy, = 0. However, because of the complexity 
the reauiting equation it was found more convenient to determine the 
| 


January, 1 


We 


- from the preceding equations, was to calculate for various values 
C, the values Per until a was reached in which 

for three | successive values of differing by 0. 1, , the value of P. cr Corres 
refinement, the relation was approximated by a second degree 


polynomial for these three points, and the final value of Poy was found aval 

the minimum of this polynomial. _ It appeared, however, that this additional _ 


_ procedure | was of more value in determining the final value of of Cc than of a 
near r the minimum does not vary appreciably. 


& 


dimensionless forms for expressing the critical 
was used by and Fung,( while 


Comparison with Available Experimental Data 


Figures 2 3 show the comparison of the with 
calculated values given in Table 1 and with other published theoretical 


a can be seen that the method used here gives a fairly good approximation 


the loads for the range X=3.5 to A =7.0. For 
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SPHERICAL SHELLS 
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5 | 0.026 63 


.0 | 0.017 47° 
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lan and Fung 


Constant “Volume 


Kaplan and Fung: 


Air Pressure 


The experimen- 
“tal results | show such a | wide scatter that an : accurate comparison is of course 
not possible. The prediction by this method in the } range of h mentioned « can - 


ae... present method, as can be noted, compares well with previously pub- . 


ae evidence | tends to show that buckling of deep spherical shell mm 
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y larger. The difference be 4 
i 
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xperimenta 
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"ESS 


4, 
essure 


Air Pr 


|) | 
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shell, clamped around the edge, lies in the different q 


is illogical to suppose, as the theory predicts, that the buckling load of a 
a shell will increase as the distance from the centre to the clamped boundary — - 
ae ce becomes larger. Hence, as an approximation, the hypothesis can be advanced ~ 
| that the buckling pressure for any depth of shell will be the same as that for 
the shallow shell, , with the same thickness and curvature, whose ey -4 
‘distance ce is such that the buckling | pressure is a minimum. 
_ It can be seen from Figure 2 that the minimum buckling ‘pressure so de- 


“fin ed gives a reasonable ee for the buckling loads for all the ex- 


cal, in case the actual 1 buckling pressure is reduced | to the dimen en- 


sionless form Per the maximum radius of curvature oft shell in 


Per, 
mental Theoretical 


e e i 0. 218, 0.457 
Be: It can be seen that the theoretical value, based on 1 the aii that 
= 0.3, agrees reasonably well with three of the four experimental results. 
discrepancy in the result was due an in the 


of the 


This paper is based on material submitted in partial | fulfillment of the 
requirements for the degree of Sc.D. at the Massachusetts Institute of Pr 
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SIMPLIFICATION OF DIMENSIONAL  ANALYSIS* 
C. | Bowman] and | Vaughn E E. Hansen 


combined | to form dimensionless parameters which can be grouped into mean - 


‘dimensionless parameters using determinantes to form simultaneous equa + 
tions is a rather tedious procedure. Considerable simplification results when — 
: the following steps are used singly o or in combination: (1) Special arrange- *- 
ment of the variables in the determinant form, (2) Selection of repeating \ vari- 
ables to represent the dynamic, kinematic and geometric variables. 
as Before using dimensional analysis the engineer must be familiar with | the 
-mechanics of the field in which he is working. He must realize dimensional | 
‘ analysis: provides a partial solution without revealing the inner phenomena 
which may occur. . It is a systematic treatment yielding a an equation which is 
dimensionally correct. The dimensionless coefficients which are formed can 


Conventional Determinant Method nod 


on _ The standard procedure of solution can be illustrated - analyzing. the edrag 


. on on a ship. _ The drag force F depends on the speed of the ship V, length of the 7 
ship L, the viscosity the mass density P, and the acceleration of gravity, 
These variables in functional as f{(F, V, L, uw, P, g) = 


= Note: Discussion open until June 1, 1959. To extend 7 closing date one month, a 
written request must be filed with the Executive Secretary, ASCE. Paper 1898 = 

a. part of the copyrighted Journal of the Engineering Mechanics Division, Proceedings ~ 
American | Society of Cir Civil ‘Engineers, vet. 85, No. 1, | 1988. 


a _be used to obtain similarity between model and prototyp 


Dept., , Spring Quarter, 1 1957, at Utah State University, Logan, Utah. 
\* Asst. Prof., Agri. Eng. Dept., , Montana State | College, Bozeman, Mont. 
-Prof., Civ. ‘and Irrig. Eng., Utah State Univ., Logan, Utah. 
& , The purpose of this paper is to show the arrangement and handling of a 
- variables rather than justifying the selection of variables. Foraclear 
treatise of the selection of variables refer to Langhaar, Dimensional oe 
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‘Seven steps are used to illustrate briefly how the is used 
+i to solve the simultaneous equations and to obtain the form of the dimension- 7 
less parameters. . After they have been formed in the aan icomaed 7 
4 the utility of the simplifying steps will be shown. awe ty at 


To obtain the dimensionless parameters, first place the variables tarsheeh | 


which is easiest to experimentally. in turn is followed by 
= easiest and so on. . The dimensions of the variables are placed on the 
side of the matrix forming the rows. The exponents of the dimensions 


‘sions. _ Positive numbers are used for exponents appearing in the semnioes 

and negative numbers for the exponents appearing in the a ‘se ae: 


- subtracting the number of dimensions from the number of variables. In the © 
ry case being considered six variables appear along the top of the matrix: and >. 
three dimensions ; along the side. Thus (n-r) = 6-3= 23. Therefore, there will 
be three dimensionless when the variables are 


a Step 3: “Obtaining Equations for the 


_ Assign K values to all variables, Ky corresponds to F, Kg corresponds to 
V, Kg corresponds to L, K4 corresponds to uw, K5 corresponds to F, and Ke 


Vv 


for on Ks and Kg in terms of Kj , Kg and K3, are 
eo established which force Kj, Kg and K3 to equal K4, K5 and Kg respectively, — 
thus | causing the F, L and T rows to be dimensionless. Since inadimen- 
a sionally homogeneous equation the dimensions of force,: length, and time 
a to Zero, following equations can be from the above 
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followin e uations are obtained: = # 
= 


ENSIONAL ANALYSIS 


+ 
= - + K6 - 0 
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Solving eq. 2, 3, and 4 for K4, Ks, Ke in terms of the 


+ 
& 


wo 


+ 
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2K - - 2/3K3, 


K5 + 1/3Kq + 


“Step 4: Setting up a Matrix for Dimensionless Parameters” at 


- Set up a new matrix for the K values listing the K values at the head of the e 
columns in their original order. The three dimensionless parameters are 7 
_ represented on the left side of the matrix by and 73. Any values 
® substituted for Ki, K2 . and K3 where K4, K5 and Kg are solved in terms of 
_ the first three K values and the parameters will be dimensionless. However, 

‘it is convenient to set Kj - = 1, Kg = 0, and Kg = 0 for 74; Ki = 0, Kp = 1 = 
Kg = 0 for 79; and Ky = 0, Kp = 0 and Kg = 1 for 73. By substituting these 
= into equations 5, 6 and 7 the values of K4, K5 and Kg for various 7 


values can be | entered in a matrix as shown below. ai So 


1 0- 1/3, vs 


0 1-920 1 

‘Forming Dimensionless | Parameters 
- The dimensionless parameters can now be read directly ly by reading from — 

the m factors to the right along the rows. The exponential factors for each of 

the variables can be read off directly and placed in the numerator or denomi-_ 
nator according to the sign. Thus obtaining the following parameters: — 

which combine into the following functional equation: 


existing are valid and equation 8 is correct, but the para- 
_ meters are not readily identifiable in their present form. The parameters _ 
should be. cleared of the fractional exponents and in the Sem: of a 


ue => (3) — ig 

® 
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‘combinations form Reynolds number, Froude number and Salen number: 


X 


which combined also form a equation 


_ To find the force on the ship, it is necessary to write eq. 9 so that the _ 


_ parameter force is explicit. 17 is the Euler number. 


This equation states that is equal 1 to a function of (Np, Np times P Vv 
The L2 can be replaced by the cross sectional area A of the boat , which is a 
determined by the design of the ship below the water line. 4 The f4(Nr, Ng) can , a 
be represented by a } REET coefficient called a drag coefficient — paneer > 


designated by substituting these values into ‘equation 5 the following 
the Matrix Method of Solution 


_ The conventional method used above is a long but useful method and can be 
used when 1 many variables are involved and the identities of the parameters — 
are unknown before hand. Due to the many steps the chance for error is e 
- great. - This method can be simplified and shortened by the arrangement | of 
the variables in step one. In this particular problem, it is realized the : a 
_ Reynolds number, Froude number and Euler number may be present in some — 
form. - Select variables which will appear only in one of the three dimension ‘ a 
“less parameters, and place these variables at the head of the first three . 
columns. The last 3 columns (this being a six-variable problem, are headed 
- by) one dynamic, one geometric, and one kinematic repeating variable mace 


all three types of variables are represented in the problem. _ The dynamic : 


_ variable represents the force dimension, the cesta represents the length, 
mi. is convenient at this point to solve for the repeating variables (V,L, P) in 4 
terms of the non-repeating variables (u, g, F) without resorting tothe K | ae 
terms which were used in the first solution. Thus — the original matrix 
form with the suggested changes in order of of listing, 
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em 4. Vennard, sony orgie Fluid | Mechanics, Third E Edition, John Wiley anc and ei be. 
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0 
‘By solving for the repeating variables in terms of the n non-repeating varia-— ad 
bles the following equations are obtained, remembering that the dimensions _ Fe 
of the variables are in the equations and not the nhumerical values 
of the variables: 


N3 | 0 2 -2 
1 From matrix the dimensionless par: param can be directly ate 


The readily recognized ratios of and Froude numbers are 


Ne = fg (Nr, 
is obtained without having to find a combination of parameters which \ will : 


a ‘This second solution proves two important points. e ar gen 


(1) The arrangement all 
the variables can simplify the solution of dimensional analysis. (2) Repeat ae 
variables must be selected which have the dimensions necessary to cancel 


a = the dimensions of the non- repeating variables. ‘When F, L and 7 appear 


kinematic, one geometric variable in the repeating variables to form the 
_ When the number of parameters have been determined by 9 (ne r) as pre-— 
“viously discussed, it also indicates there will be (n- r) variables which will - ai 
appear alone in each of the dimensionless (n- 


4 
— |. 
Using these equations to form a second matrix a e 
a for identi 
| 


ables which cannot be identified as ; part of a known number and these in poems 4 
are followed by one dynamic, one kinematic, and one geometric repeating i aN 


5.5 The parameter probably, , although ‘not necessarily, v will = 
identifiable numbers but they will be valid dimensionless eae which 
7 _ are characteristic of the problem being solved. d. 


Simplified | Dimensional Analysis 
The above observation ‘introduces a ‘new to the procedure of 
forming dimensionless parameters. By using | the original problem of the pan 
2 drag on a ship which is a function of (F, Vv, te HM, P, g) the simplification can 
_ be made quite evident. There are six variables present which are made up of | 
three dimensions, F, Li and T. Since (n-r) equals 3, there are three possible 
dimensionless parameters fo1 formed from the above variables. By 
and knowledge of the problem it is known that uw will appear in the Reynolds — 
number ana g will appear in the Froude number, thus leaving one nla 
unknown. . By selecting one dynamic, one kinematic and one geometric repeat ez 
- ing variables which are known to appear in both of the known parameters, the _ 
_ third non-repeating variable can be determined. It is known that V, L appears 
in both the Reynold and Froude numbers. "Therefore, it is safe to use them veal 


two of the required three repeating variables. 


_ The variable V represents the kinematic variables involving t time and Lis is. 


‘representative of the geometric variables involving length. A third d dynamic — 


_ The Reynold number gives the third repeating variable because this num=— : 
‘ing g. Density P which appears” in the e number must be the other ‘repeating ee 
variable. _ Thereby, the three repeating variables (V, L, ?) have been selected 
leaving the variables (u, g, F) as the non repeating variables. 
_ To form the parameter, first select one of the non- -repeating variables in 
quotient form such as yw, then combine this variable with any or allofthe | 
_ three repeating variables in such a way that the combination has zero dimen- 
sions. The following shows how thisisdone. 
we Since uw has a dimension of a divide by the acme repeating variable 
te cancel the force dimension 
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5. ‘if If only dynamic kinematic variables only on one @ dynamic ana 
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zero dimensions = N, (Reynolds number) 


Repeating the above steps for the remaining non- ~repeating vartenies 


Dividing t by eliminates force 


= zero dimensions 


—- = zero dimensions © 


2 = N, (eaters 


_ The variables are now formed into @menatenines parameters which can 
be combined to give ed 9, from which eq. i can be obtained as was done i in 


‘The first method derives the dimensionless in seven steps. 
g many problems this procedure is desirable when the user is not inti- = _ 


mately familiar with the field in which he is working. ‘The second method 

eliminates two steps from the first solution and derives the same parameters 

< with greater ease. ‘The ease by which the solution is made is based on the 


= ‘This arrangement is s dependent again on the Inowledge of the user. ah ; 
- _ The third method is very simple and can be used merely by careful selec-_ 
¥ of the repeating variables and then forcing the remaining variables into 
parameters. The parameters are identical to those found by 
the other two methods. 
_ ‘The three methods « as presented are valid for finding | dimensionless 
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‘DY NAMIC ELASTO-PLASTIC OF RIGID 7 


_ 


_ Discussion by A. T. Matthews 


wk to the dynamic analysis of structures. ee 


~ To illustrate this method of expansion into normal vibrational modes the 
- dynamic response of a two-hinged rigid frame was obtained and carried up | 
wf to the start of the second elasto-plastic phase. Regarding the assumptions | 
and certain restrictions which are implied in this } analysis th the ee ob- 
‘servations may be made: 


a In the example ¢ given, the location of the plastic neee is assumed to be 
known, and it coincides with the position of maximum moment under static 
— loading. ; a is true that hinges will form where the stress is maximum if the a 
- cross-section of the members is constant, and this in turn implies that hinges , 
form at the location of maximum moments. . However, under dynamic loading 
the effect of inertial forces may be such that maximum moments may occur 
_ at positions different from the static case. It is conceivable for certain types” 


_ of two hinged frames that the maximum dynamic loading may not occur at the 7 7 
_ knee, even under the loading assumed by the author. However, using the a - 
author’ s method it is possible to determine the actual position of the maxi aan 


mum dynamic moments at the end of the first elastic phase. _ Therefore the -. - . 


the effect of the axial forces in addition to moments in ee eran eal 

2. ‘Subsequent elasto- -plastic phases 


lose Under transient loads of sufficiently long duration repented elastic and 
- elasto-plastic sequences may be anticipated, with hinge formation phenomena 
be, as discussed above. _ ‘The author’ ’s method requires a new set of modes when- 
“ever a new plastic hinge occurs. In general the computational work of deter-— 
mining responses becomes impenction: when many hinges form and re- -form 
-—- requires knowledge of moments ey points at all times during the ritsttigaeal 
_ and determination of these moments in the elasto-plastic phase requires _ iw! 
summation of the second derivatives of displacements in each plastic mode. — 
Poor convergence of such a series, as acknowledged by the author, ae 
a. Proc. Paper 1693, July 1958, by Frank L. DiMaggio. “fans : 
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‘the “moment determination the first -plaatic phase of the motion. 

Therefore, this method is applicable to cases where the duration of the load 7 

short compared to the period of the | lowest elasto- “plastic mo mode and 
practical engineering results for such cases. 

ie In a recent paper by H. H. Bleich,(1) an alternative method is ties in 

which the elastic modes: are > retained throughout the entire time his- 
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THE RAPID INTORSION 
THE RAPID DESIGN OF BEAMSIN TORSION ae 
"Cedric Marsh? 


‘The design of beams in which occurs, usually due 13 
of the applied load, is either crudely approximated or, if studied accurately, _ 7 
es the use of hyperbolic functions. _ This paper | shows how simple — 


si 
a can be obtained. 


usual to employ shapes of the closed type, such as tubes or box beams, in ‘a 
which case analysis is relatively direct. the twist arises due to the 
7 4 centricity of lateral loads, whose p primary effect is bending, the shapes used 
oa will normally be of the ‘standard I form, possibly augmented by channels fix ll 
to one, or both flanges. This paper deals particularly with such shapes. ae 
A structural shape resists twisting forces by a combination of torsional _ 
‘Tigidity, defined by the term GJ, and torsion-bending rigidity defined by the — 
‘ aa term EH. In these terms, Jis the torsion constant, H is the torsion-bending — 
In an nt shape, that portion of the torque resisted the torsional rigidity 
_ creates shear stresses as shown in Fig. 1a, while that portion resisted by 
torsion-bending creates shear stresses in the flanges (Fig. 1b) which are ac- 4 
companied longitudinal stresses as illustrated in Fig. 1c. These longitudi-— 
nal stresses, as they are in addition to the longitudinal stresses caused by 7 4 
a simple bending, are of first importance, and the following approximate 7 ae 
sis concerns itself specifically with these stresses. This treatment may also A ’ 
be used to give the deformation. Only in rare cases will the shear stress due 
to pure torsion in open sections govern before deformation becomes unaccepta-— = 


Note: Discussion open until June i, , 1959. To extend the closing date one month, a = * - 

written request must be filed with the Executive Secretary, ASCE. - Paper 1923 is 4 

part of the copyrighted Journal of the Engineering Mechanics Division, Proceedings ; 
f the American jan Society of Civil Vol. 85, No. EM 1, 1969. 
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IN TORSION 


the b of: a in iat terms familiar to engineers 
it is convenient to deal —_ torsion- bending and torsional rigidities inde- __ 


‘There is an exact between the behaviour of a beam in bending and 
i a beam resisting torques by torsion- ~bending © only, as it does in in stocky | beams. 


_ If the maximum bending moment and deflection in a aa are _ by the 


where C n and are the appropriate constants, 


W is the total load, 


L is the span, 


where his is the distance from the centroid to the extreme fibre, then, for a _ 
subjected to a total torque, ia the same manner as W, the 
maximum torque-moment and rotation are: 


EH 
and the maximum longitudinal stress eee by this torque is is: 


were. isa of the ‘shape, ‘called here “ “unit warping,” 


In this case | = ‘1/8 and m, 
built in the e carrying a torque, - 


— 
| 
A 
— 
~, 
— 
q iim 


_ if the load W is applied such that the distance from the shear center to the 
ve Hine of action of the load is e, the torque created on the member is We. Thus, 7 
_ it follows that, neglecting torsional rigidity, the ratio between the maximum > 


stress due to ) the torque, We, and the stress caused by direct bending due ‘all 


will be called 


If, as in 1 slender beams, the torsion- a can be and 

_ only the torsional rigidity is considered, the relevant laws governing stress _ 
7 _ and rotation will be analogous” to o those governing shear stress & shear deflec- 


a a aie For example, a load W applied at the quarter point of a simply supported ; 


by re GV “represents the 2 shear rigidity c of the ie beam. A to torque T applied at 
quarter point of a Simple b beam gives end of T and T anda ro- 


tation of — shear deflections are not usually calculated, 


the formulas may be unfamiliar, the equations are simple enough to be com- 


a “puted directly. (Appendix b) isa ae of some typical expressions for these 
deflections), 
Combined Torsion and 


_ beam creates reactions of Ry Mammy W, and a shear deflection of = 


approximately by proportioning them in the ratio of the 


‘This means that uf, under a a total spplied torque T , the beam 


— 

— 
— 

| 

a 
| 

= an : Considered independently it is seen that the two modes of resisting torque a’ 
te follow relatively simple laws. The interaction of torsional andtorsion- 

bending rigidity, however, involves a differential equation whose solutionis 

in hyperbolic functions. To avoid this, it is here proposed that the partition- | 

| 

64 = Ct , iftorsion-bending rigidity is neglected, | 

4 _ the ratio of the torque carried by torsion-bending, Ty, to that carried by 


BEAMS 
"pure torsion is 


acu? = a, 


“a The longitudinal stress due the wai T is then, , from (1), (2) and (4): 


und the ratio between the longitudinal stress due to twisting c caused by a 


eccentrically, applied load and that due to bending, is: Omi 


Se 
bei + Ca) 


is: 


maximum rotation caused wy its e, is: 


and the total movement of the load i is: 


a In shapes symmetrical about both axes, , such as I beams, the shear center 
7" In channel mages the shear center, lies on the x-x axis ata distance e- 
behind we, the median line of the (Fi 


flanges, xX is the distance from the ‘me@in line of the web to the centroid, and — 
i is the radius of gyration about the Xx axis. (See ref. . 1 for the other shear 

Va When two shapes are combined, the shear center lies at a point on th the axis } 


joining | the shear centers ‘of f the individual ‘Shapes: 3) such that: 


ats 


Let C 5 / C, = C, an e,* 
a 
+ 
| 
4 
mation noe of th 
fhe donth Of the chang red hotwoon the mas 
| — 
f tia of the individual shape at 
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IN TORSION 
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Torsion Cuan 


For a thin wall of breadth | b thickness the to constant is: 


_ The sum of these terms for all the individual walls will give a conservativ ve 
approximation to the torsion constant of a a shape. A more exact figure takes” -_ 
account of the influence of the roots fillets. For many standard shapes the _ 

values of the torsion sion constants 3) and the calculation is is not 


_hecessary. 


-Bending Constant 
torsion- -bending constant of an I-shape is: 
= 


is the m moment of inertia about the y-y axis, : 
d is the depth, measured between the median lines o of the 
Fora channel ‘(he 2)= Nay. 


sh where ry is the radius of gyration about the y y- -y axis. is. i 
two ‘Shapes are combined the torsion- bending constant is equal to the 

of the individual torsion- bending plus the term: 


ie 


- for each individual shape, where ‘ach is the distance from the shear center | of 
_ the individual shape to that of the complete section, and I is the moment of a 


section, ‘which would have plane during simple to “warp” 4 
out of plane under the action of a torque. The term “unit warping,” D, is oe 
4 if such warping is s opposed, ‘it is s related to the 1 maximum n stress created. iil - a 
Ett The calculation of D will be illustrated bd by the example of an I- ~shape with a a 
_ channel f fixed to each flange (Fig. 4a) 
‘The part ABC of the flange in Fig. 4a is ; opened ‘out (Fig. 4b). This line is 
- then treated as a beam subjected to a uniform load equal to the Bin tnd ; 
_ distance from the shear center of the - complete section to the median line =. 4 ; 
the wall in question (Fig. 4c). The total load on the beam is then reacted at Som 
} the end representing the free edge of the flange by a force D. . The shear stress 
_ in this fictitious beam (Fig. 4d) represents the distribution of stress in the a 
flange due to torsion-bending, the maximum stress occuring at the free edge. 
ih the other of the flange the stress distribution of 


| 
(10) 
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the products of the breadth of a wall and its distance from the shear center, ea 


such sum comprising only the material in one half of a flange assembly. Ree 


Example | (Fig. 5). 


x simply aluminum ‘comprised of a x4. 33" x a? 
- If. I, with a 6"! x 1.92'' x 2.83 lb/ft. channel fixed to the top flange, car- 
7 ries a cental load W applied 2 2'' above the back of the channel web, over at 
span of Find the maximum (For aluminum E= 


q 107 psi, G = 4x 108 psi). 
q 


that this term is of is the sum net 


9 


= 6. -38 in 


shear center is at the centroid. 
wary warping ‘constant is: > 


=9-2x0.29 


2 in. 


The shear en lies at a distance above the median line of the channel 


-C-Shape e 


» 


= 6- 0.4 = 5.6 in., 


= 0. 41 (5.6/2 x 2. 34)? = = 0.59 in, 

The torsion- -bending 


= 


For the: 4 
| 


~ 


0.51. -0. 1 = 0.41 
al 


1 


0.41) 


= 


0. 
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BEAMS IN TORSION 


in 


§ 


on 


» 
(| 


The distance from of the I-shape to the shear center of the 


el = 2 


aiid 


1938/7. 6 in in.4 
18.21 ind 


0. 548 in.4 


“For a simply beam with a centr 


=1/4 
Hf 


Ca) 
a +Ca) =e x 7. 3) 
7 
yer stress due to direct bending is: 
Wx80 


B54 x 1.62 


‘giving a total stress of: 


— 
— 
— ‘ 
| 


at the center is: 7 


Wx 3x 803 


48x 1.62x 10” x 193 


e 59 + + 3. x 1.02) x 10~ 
= 8. 65 x in. 


‘The stress calculated by this meteel is — ‘accurate for stocky beams. 
For a unigneen load of constant eccentricity the error is negligible for all _ 
spans. Fora , CORES load, the error is small for values of “a” up to 1 _ 
: = -10% for “ a” =2 , (as ap percentage of the total stress it is appreciably — 
‘than this). As a” becomes larger, deformation is more likely to govern 
sign and errors in the computed stress are less significant. 
Rotations computed in the above way are accurate for all values of “a” 7 
the torque is distributed. A maximum error of approximately +10% 
- occurs for a central torque when “ a” lies between 3 and : 4, , the error approach- 
ing zero as “a ” approaches zero or infinity. 
demonstrate the degree of accuracy, c) is an example taken 


‘pared with the. exact answers. 


As the values of torsion constants are published, and the torsion- “bending 
oo of shapes can be computed from other published properties, the = 
4 analysis of beams in torsion, using existing engineering formulas, is seen to | 


APPENDIX. 


a) The fundamental equations governing rotation in a member, due to ap- 
plied torques opposed wholly torsion - bending are ‘given below 


| 


48 (1+ Ca)EH ~ x W radians. 
4. 
4 — 
| 
[| 
— — 
— 
a 
— 
ss Where t is the app o the applied — 
where I is the va i> 


“EI EH 
is here c: called the -moment, analogous tot the 


where > i is the analogous to the deflection, 
Equations gov governing rotation created by pure 
“torsional rigidity, with those due to shear, are given 


“below. 


ret 
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This is one given in Stresses in 
ae Simply supported 14WF74 spanning 20 ft. with a central vertical load of 
_ Effective d d = 13.4 in, 
Flange breadth b = 10.07 in.. d a b/4 = 13.4 x 10. 01/4 = 33.1 7 in. 
= 133.5 in. 42 000. in. 
199.525 = 6,000 
x 112.3 in., U=H/D = 6000/33. 


J 924 in.4 


(2 


(sca) = 


q 
The max rigidity, are given for some 
| 
— | 


(112.3 x2 


_ Maximum Stress = 
18500 ps psi 
f. 19200 psi by exact method 
x 2403 


Maximum Rotation = 8 Ca) EH a8 * 
031 radian 
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constant in formula for shear 


“unit warping,” : a a to give 
the depth of an I or channel measured between the ‘median 


lines of the flanges, in. > 


the e elastic psi., 


torsion - stress, 


w. 
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constant for a shape, in.* 
distance from the centroid to ) the extreme fibre of a shape, | in. aaa 

= moment of inertia of a shape, in. ; 


= torsion constant for a shape, in.’ 


span of a beam, in., 7 
Ag 


maximum moment t, Ib. in, 
radius of gyration, 
= ry modulus ofa shape, in.3 
te = total applied t torque, |b. i 
_= torque carried by | torsional rigidity, Ib. in. 


= torque carried by torsion- -bending rigidi » Ib. 


shear constant for a 


q 


maximum 

rotation, radians. 


w 


— 

— 
— 

ir 
— 
ty 

= a distance measured from the centroid, in., ie 


